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The effects of the propagation of particles which have a finite Hfe-time and an 
according width in their mass spectrum are discussed in the context of transport 
descriptions. In the first part the coupling of soft photon modes to a source of 
charged particles is studied in a classical model which can be solved completely 
in analytical terms. The solution corresponds to a re-summation of certain field 
theory diagrams. The general properties of broad resonances in dense finite tem- 
perature systems are discussed at the example of the p-meson in hadronic matter. 
The second part addresses the problem of transport descriptions which also ac- 
count for the damping width of the particles. The Kadanoff-Baym equation after 
gradient approximation together with the <I>-derivable method of Baym provides a 
self-consistent and conserving scheme. Memory effects appearing in collision term 
diagrams of higher order are discussed. We derive a generalized expression for 
the nonequilibrium kinetic entropy flow, which includes corrections from fluctua- 
tions and mass-width effects. In special cases an J?-theorem is proved. Memory 
effects in collision terms provide contributions to the kinetic entropy flow that in 
the Fermi-liquid case recover the famous bosonic type In T correction to the 
specific heat of liquid Helium-3. 



With the aim to describe the collision of two nuclei at intermediate or 
even high energies one is confronted with the fact that the dynamics has to 
include particles like the delta or rho-meson resonances with life-times of less 
than 2 fm/c or equivalently with damping widths above 100 MeV. Also the 
collisional damping rates deduced from presently used transport codes are of 
the same order, whereas typical mean temperatures range between 50 to 150 
MeV depending on beam energy. Thus, the damping width of most of the 
constituents in the system can by no means be treated as a perturbation. As 
a consequence the mass spectrum of the particles in the dense matter is no 
longer a sharp delta function but rather acquires a width due to collisions 
and decays. One thus comes to a picture which unifies resonances which have 
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already a width in vacuum due to decay modes with the "states" of particles in 
dense matter, which obtain a width due to collisions (coUisional broadening). 

The theoretical concepts for a proper many body description in terms oLa 
real time nonequilibpium field theory have already been devised by SchwingerEI, 
KadanofF and BaymH, and KeldyshO in the early sixties, extensions to relativis- 
tic plasmas followed by Bezzerides and DuBoisQ. First investigations of the 
quantum effects on the Boltzmann collision term were given by DanielewiczEl, 
the principal conceptual problems on the level of quantum field theory were 
investigated by Landsmann El, while applications which seriously include the 
finite width of the partirlf't; in tp??^n"r^^ descriptions were carried out only in 
recent times, e.g. refs. dffliJ'E^IOllj'Ej. For resonances, e.g. the delta reso- 
nance, it was natural to consider broad mass distributions and ad hoc recipes 
have been invented to include this in transport simulation models. However, 
many of these recipes are not correct as. they violate some basic principles like 
detailed balance (see discussion iiL_ref . El) , and the description of resonances in 
dense matter has to be improvedllj. 

In this contribution the consequences of the propagation of particles with 
short life-times are discussed. First we address a genuine problem related 
to the occurrence of broad damping width, i.e. the soft mode problem. At 
the classical level we investigate the coupling of a coherent classical field, the 
Maxwell field, to a stochastic source described by the Brownian motion of a 
charged particle. In this case the classical current-current correlation func- 
tion, can be obtained in closed analytical terms and discussed as a function 
of the macroscopic transport properties, the friction and diffusion coefficient 
of the Brownian particle. The result corresponds to a partial re-summation of 
photon self-energy diagrams in the real-time formulation of field theory. Sub- 
sequently the properties of particles with broad damping width is illustrated at 
the example of the p-meson in dense matter at finite temperature. In the final 
part we discuss how particles with such broad mass-width can be described 
consistently within a transport theoretical picture. 

We are going to argue that the Kadanoff-Baym equations in the fisst gradi- 
ent approximation together with the ^-functional method of Baymllil provide 
a proper frame for kinetic description of systems of particles with a broad 
mass-width. To this end, we discuss relevant problems concerning charge and 
energy-momentum conservation, thermodynamic consistency, memory effects 
in the collision term and the growth of entropy in specific cases. For simplicity 
we concentrate on systems of non-relativistic particles. Generalization to sys- 
tems of relativistic particles and hosonic mean fields can be straight forwardly 
done along the lines given in ref. c3. 
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Figure 1: Left: Current-current correlation function in units of < t;^ > as a function of 
time (in units of l/Tx) for different values of the photon momentum = 3k^T^/< > 
with k = 0,1,2,3. Right: Rate of real photons dfi N/ (duidt) in units of 4776^ ^ii^^ /3 for 
a non-relativistic source for Fx =50,100,150 MeV; for comparison the IQF results (dashed 
lines) are also shown. 



1 Bremsstrahlung from Classical Sources 

For a clarification of the soft mode problem we discus an example in classical 
electrodynamics. We consider a stochastic source, the hard matter, where the 
motion of a single charge is described by a diffusion process in terms of a 
Fokker-Planck equation for the probability distribution / of position x and 
velocity v 

|/(->-^^) = (^r^|^ + r4.-.^) (1) 

Fluctuations also evolvain time according to this equation, or equivalently by 
a random walk process E3, and this way determine correlations. This charge is 
coupled to the Maxwell field. On the assumption of a non-relativistic source, 
this case does not suffer from standard pathologies encountered in hard thermal 
loop (HTL) problems of QCD, namely the coUinear singularities, where vq w 1, 
and from diverging Bosc- factors. The advantage of this Abelian example is 
that damping can be fully included without violating current conservation and 
gauge invariance. This problem is related to the Landau-Pommeranchuk- 
Migdal effect of bremsstrahlung in high-energy scattering 

The two macroscopic parameters, the spatial diffusion D and friction F^; 
coefficients determine the relaxation rates of velocities. In the equilibrium limit 
(t oo) the distribution attains a Maxwell-Boltzmann velocity distribution 
with the temperature T = mlv"^^ /3 — mDTx- The correlation function can 
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be obtained in closed form and one can discuss the resulting time correlations 
of the current at various values of the spatial photon momentum q, Fig. (de- 
tails are given in ref.Ej). For the transverse part of the correlation tensor this 
correlation decays exponentially as ~ e~'""='^ at q = 0, and its width further de- 
creases with increasing momentum q — \q\. The in- medium production rate is 
given by the time Fourier transform r — > w, Fig. (right part). The hard part 
of the spectrum behaves as intuitively expected, namely, it is proportional to 
the microscopic collision rate expressed through (cf. below) and thus can 
be treated pertubatively by incoherent quasi-free (IQF) scattering prescrip- 
tions. However, independently of the rate saturates at a value of ~ 1/2 
in these units around lo T^, and the soft part shows the inverse behavior. 
That is, with increasing collision rate the production rate is more and more 
suppressed! This is in line with the picture, where photons cannot resolve the 
individual collisions any more. Since the soft part of the spectrum behaves 
like uj/Tx, it shows a genuine non-perturbative feature which cannot be ob- 
tained by any power series in T^. For comparison: the dashed lines show the 
corresponding IQF yields, which agree with the correct rates for the hard part 
while completely fail and diverge towards the soft end of the spectrum. For 
non-relativistic sources (^v^) <C 1 one can ignore the additional g-dependence 
(dipole approximation; cf. Fig. ^ and the entire spectrum is determined by 
one macroscopic scale, the relaxation rate F^. This scale provides a quenching 
factor 

by which the IQF results have to be corrected in order to account for the finite 
collision time effects in dense matter. 



The diffusion result represents a re-summation of the microscopic Langevin 
multiple collision picture and altogether only macroscopic scales are relevant 
for the form of the spectrum and not the details of the microscopic collisions. 
Note also that the classicaLijeault fulfill the classical version (h 0) of the 
sum rules discussed in refs. EZrEj. 



2 Radiation on the Quantum level 

We have seen that at the classical level the problem of radiation from dense 
matter can be solved quite naturally and completely at least for simple ex- 
amples, and Figs. ^ display the main physics. They show, that the damping 
of the particles due to scattering is an important feature, which in particular 
has to be included right from the onset. This does not only assure results 
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that no longer diverge, but also provides a systematic and convergent scheme. 

On the quantum level such problems require techniques beyond the standard 
repertoire of perturbation theory or the quasi-particle approximation. 

In terms of nonequilibrium diagrammatic technique in Keldysh 
notation, the production or absorption rates are given by pho- 
ton self-energy diagrams of the type to the right with an in- 
and out-going photon line (dashed). The hatched loop area 
denotes all strong interactions of the source. The latter give rise to a whole 
scries of diagrams. As mentioned, for the particles of the source, e.g. the nuclc- 
ons, one has to re-sum Dyson's equation with the corresponding full complex 
self-energy in order to determine the full Green's functions in dense matter. 
Once one has these Green's functions together with the interaction vertices at 
hand, one could in principle calculate the required diagrams. However, both 
the computational effort to calculate a single diagram and the number of di- 
agrams increase dramatically with the loop order of the diagrams, such that 
in practice only lowest-order loop diagrams can be considered in the quantum 
case. In certain limits some diagrams drop out. We could show that in the 
classical limit, which in this case implies the hierarchy cj, |q|, F <C T ^ to to- 
gether with low phase-space occupations for the source, i.e. f{x,p) <C 1, only 
the following set of diagrams survives 

■O- ^ ■'<$>'■■ ^ - <tf 't> ' - 

In these diagrams the bold lines denote the full miclcon Green's functions 
which also include the damping width, the black blocks represent the effective 
nucleon-nucleon interaction in matter, and the full dots the coupling vertex 
to the photon. Each of these diagrams with n interaction loop insertions just 
corresponds to the n*'* term in the corresponding classical Langevin process, 
where hard scatterings occur at random with a constant mean collision rate 
r. These scatterings consecutively change the velocity of a point charge from 

Do to to t;2, In between scatterings the charge moves freely. For such 

a multiple collision process the space integrated current- current correlation 
function takes a simple Poisson form 

in''''-+oc J d'xid'x2{r{xi,t - ^)j^^{x2,t+ ^)) 

°° IFtI" 
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with V = {l,v). Here (. . .) denotes the average over the discrete colhsion se- 
quence. This form, which one writes down intuitively, agrees with the analytic 
result of the quantum correlation diagrams in the limit n <IC 1 and T <^T. 
Fourier transformed it determines the spectrum in completely regular terms 
(void of any infra-red singularities), where each term describes the interference 
of the photon being emitted at a certain time or n collisions later. In special 
cases where velocity fluctuations are degraded by a constant fraction a in each 
collision, such that (dq • Vn) = a" {vq ■ Vq), one can re-sum the whole series in 
Eq. (Q) and thus recover the relaxation result with 2Tx {v'^)— T {{vq — f i)^) 
at least for q — and the corresponding quenching factor (g). Thus the clas- 
sical multiple collision example provides a-nuite intuitive picture about such 
diagrams. Further details are given in ref.113. 

The above example shows that we have to deal with particle transport 
that explicitly takes account of the particle mass-width in order to properly 
describe soft radiation from the system. 

3 The p-meson in dense matter 

Another example we like to discuss concerns properties of the p-meson and 
their consequences for the p-decay into di-leptons are discussed. In terms of 
the nonequilibrium diagrammatic technique, the exact production rate of di- 
leptons is given by the following formula 




= fp{m,p,x,t) Ap{m,p,x,t)2mT'"' (m). (5) 

Here F'' (to) oc is the mass-dependent electromagnetic decay rate of 

the p into the di-lepton pair of invariant mass to. The phase-space distribution 
fp{m,p,x,t) and the spectral function Ap{m,p,x,t) define the properties of 
the p-meson at space-time point x, t. Both quantities are in principle to be de- 
termined dynamically by an appropriate transport model. However till to-date 
the spectral functions are not treated dynamically in most of the present trans- 
port models. Rather one employs on-shell (5-functions for all stable particles 
and spectral functions fixed to the vacuum shape for resonances. 

As an illustration, the model case is discussed, where the p-meson just 
strongly couples to two channels, i.e. the t:^tt~ and ttTV <-^- pN channels, 
the latter being relevant at finite nuclear densities. The latter component is 
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representative for all channels contributing to the so-called direct p in trans- 
port codes. For a first orientation the equilibrium properties^ are discussed in 
simple analytical terms with the aim to discuss the consequences for the im- 
plementation of such resonance processes into dynamical transport simulation 
codes. 

Both considered processes add to the total width of the p-meson 

Ttot{in,p) ^Tp^^+^-{m,p) + Tp^^NN-^{m,p), (6) 

and the equilibrium spectral function then results from the cuts of the two 
diagrams 




Ap{m,p) = / + i \=^ . (7) 



m2-m2-ReS^)^ + m2r2 



tot 



In principle, both diagrams have to be calculated in terms of fully self-consistent 
propagators, i.e. with corresponding widths for all particles involved. This 
formidable task has not been done yet. Using micro-reversibility and the prop- 
erties of thermal distributions, the two terms in Eq. contributing to the di- 
lepton yield (||) , can indeed approximately be reformulated as the thermal aver- 
age of a TT+TT" — > p ^ e"'"e~-annihilation process and a TriV pN e^e~ N- 
scattering process, i.e. 



dmdi 

+ hfN v^N <j{7rN -^pN^ e+e-N))^ , (8) 

where and /at are corresponding particle occupations and Utttt and v^^n sue 
relative velocities. However, the important fact to be noticed is that in order to 
preserve unitarity the corresponding cross sections are no longer the free ones, 
as given by the vacuum decay width in the denominator, but rather involve 
the medium dependent total width (|^). This illustrates in simple terms that 
rates of broad resonances can no longer simply be added in a perturbative way. 



"Far more sophisticated and in ppstp-tmljMfpi-jonsistcnt equilibrium calculations have already 
been presented in the literature P^ r^rr^ ' PHr^ . It is not the point to compete with them at 
this place. 
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Di-lepton rates from thermal p- mesons (T = 110 MeV) 




m^r 2m^ nip 1 GeV 2m^ rUp 1 GeV 



Figure 2: e+e~ rates (arb. units) as a function of the invariant pair mass m at T = 110 
MeV from tt+tt" annihilation (dotted hne) and direct p-meson contribution (dashed line), 
the full line gives the sum of both contributions. Left part: usinethe free cross section recipe, 
i.e. with Ftot = ^+^- ; right part: the correct partial rates (W). Ap is given by the thick 
line. The calculations are done with Fp«^^(mp) = 150 MeV and T^^^j^jy-i (rrip) = 70 
MeV. 



Since it concerns a coupled channel problem, there is a cross talk between the 
different channels to the extent that the common resonance propagator attains 
the total width arising from all partial widths feeding and depopulating the 
resonance. While a perturbative treatment with free cross sections in Eq. (||) 
would enhance the yield at resonance mass, m — nip, if a channel is added, cf. 
Fig. 2 left part, the correct treatment even inverts the trend and indeed 
depletes the yield at resonance mass, right part in Fig. 2. Furthermore, one 
sees that only the total yield involves the spectral function, while any partial 
cross section refers to that partial term with the corresponding partial width in 
the numerator! Unfortunately so far all these facts have been ignored or even 
overlooked in the present transport treatment of broad resonances. Compared 
to the spectral function both thermal components in Fig. 2 show a significant 
enhancement on the low mass side and a strong depletion at high masses due 
to the thermal weight / oc exp(— po/T) in the rate (^). 

As an example we show an exploratory study of the interacting system of 
TT, p and oi-mesons described by the ^-functional 
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Dilepton rates (arbtr. units), T=150 MeV 



rho-meson spectral function, T=150 MeV 
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Figure 3: left part: contour plot of the self-consistent spectral function of the p-meson as a 
function of energy and spatial momentum; right part; thermal di-lepton rate as a function 
of invariant mass at p = 300 MeV / c 

TT TT TT 



(9) 



(cf. section ^ below). The couplings and masses are chosen as to reproduce 
the known vacuum properties of the p and ai meson with nominal masses 
and widths vip = 770 MeV, m^i = 1200 MeV, Vp = 150 MeV, Ta^ = 400 
MeV. The results of a finite temperature calculation at T = 150 MeV witJi 
all self-energy loops resulting from the ^-functional of Eq. computed 
with self-consistent broad width Green's functions are displayed in Fig. 3 (cor- 
rections to the real part of the self-energies were not yet included). The last 
diagram of <& with the four pion self-coupling has been added in order to supply 
pion with broad mass- width as they would result from the coupling of pions to 
nucleons and the A resonance in nuclear matter environment. As compared to 
first-order one-loop results which drop to zero below the 2-pion threshold at 
280 MeV, the self-consistent results essentially add in strength at the low-mass 
side of the di-lepton spectrum. 

4 Quantum Kinetic Equation 

The two above-presented examples unambiguously show that for consistent 
dynamical treatment of nonequilibrium evolution of soft radiation and broad 
resonances we need a transport theory that takes due account of mass-widths 
of constituent particles. A proper frame for such a transport is provided by 
Kadanoff-Baym equations. We consider the Kadanoff-Baym equations in the 



9 



first-order gradient approximation, assuming that time-space evolution of a 
system is smooth enough to justify this approximation. 

First of all, it is helpful to avoid all the imaginary factors inherent in 
the standard Green function formulation (G*^ with i,j G { — h}) and intro- 
duce quantities which arc real and. in the quasi-homogeneous limit, positive 
and therefore have a straightforward physical interpretation, much like for the 
Boltzmann equation. In the Wigner representation we define 



for the generalized Wigner functions F and F with the corresponding four- 
phase-space distribution functions f{X,p) and Fermi/Bose factors [l^f{X,p)], 
with the spectral function A{X,p) and the retarded propagator . Here and 
below the upper sign corresponds to fermions and the lower one, to bosons. 
According to relations between Green functions G'-' only two independent real 
functions of all the G^^ are required for a complete description. Likewise the 
reduced gain and loss rates of the collision integral and the damping rate are 
defined as 



where S*-' are contour components of the self-energy, and is the retarded 
self-energy. 

In terms of this notation and within the first-order gradient approximation, 
the Kadanoff-Baym equations for F and F (which result from differences of 
the corresponding Dyson's equations with their adjoint ones) take the kinetic 
form 



F{X,p) = A{X,p)f{X,p) = (T)iG-+(X,p), 
F{X,p) = A{X,p)[l T f{X,p)] = iG+-(X,p), 
A{X,p) = -2\mG^{X,p) = F±F = i (G+- - G-+) 



(10) 
(11) 



r,^{x,p) = (T)is-+(x,p), ro„t(^,p) = is+-(x,p), 

TiX,p) = -2ImS^(X,p) = rout(^,p) ±ri„(X,p), 



(12) 
(13) 



PF- {ri„,ReG^} = G, 
PF-{rout,ReG^} =tC 



(14) 
(15) 



with the drift operator and collision term respectively 





C{X,p) = Ti4X,p)F{X,p) - Tont{X,p)F{X,p). 



(17) 
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Within the same approximation level there are two alternative equations 
for F and F 

MF - ReG^Fin = J ({F, F} - {T,^, A}) , (18) 

MF-ReG^Fout = J ({r,-^} - {rout, A}) (19) 

with the "mass" function M{X,p) = po - /2m - ReE^(X,p). These two 
equations result from sums of the corresponding Dyson's equations with their 
adjoint ones. Eqs. (|l|) and (|I^ can be called the mass-shell equations, since 
in the quasiparticle limit they provide the on-mass-shcU condition M = 0. 
Appropriate combinations of the two sets (p^)-(p^) and (p^)-(p^ provide Jia 
with retarded Green's function equations, which are simultaneously solved ErQ 
by 

f A- ^ 

CR - \ ^ J M'^ + r V4 ' ^20^ 

- M{X,p)+inX,p)l2^\^^Cn_ M 

I ~M2 + r2/4- 

With the solution (|2^) for G^ equations (|l^ and (|l^) become identical 



to each other, as well as Eqs. ([L8|) and ([191). However, Eqs. (|14D-(15) still 
are not identical to Eqs. (^8|)-(^9^, while they were identical before the gra- 
dient expansion. Indeed, one can showN that Eqs. (ll|)-(|ll) differ from Eqs. 



(|18|)-(|19D in second order gradient terms. This is acceptable within the gradi- 
ent approximation, however, the still remaining difference in the second-order 
terms is inconvenient from the practical point of view. Following Botermans 
and MalflictI, we express Fin = E/ -(- 0{dx) and Fin = F(l T /) + 0{dx) 
from the l.h.s. of mass-shell Eqs. ([T8| ) and (|T9| ), substitute them into the 



Poisson bracketed terms of Eqs. (14) and (p_q), and neglect all the result- 



ing second-order gradient terms. The so obtained quantum four-phase-space 
kinetic equations for F — f A and = (1 =p f)A then read 

P(/A)-{F/,RcG«} =G, (21) 
P((lT/)^)-{r(lT/),RcG«} =tC. (22) 

These quantum four-phase-space kinetic equations, which are identical to each 
other in view of retarded relation (^0[), are at the same time completely iden- 
tical to the correspondingly substituted mass- shjeU Eqs. (H) and (|l|). 

The validity of the gradient approximation^ relies on the overall smallness 
of the collision term C = {gain — loss} rather than on the smallness of the 
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damping width T. Indeed, while fluctuations and correlations are governed by 
time scales given by F, the Kadanoff-Baym equations describe the behavior of 
the ensemble mean of the occupation in phase-space F{X^p). It implies that 
F{X,p) varies on space-time scales determined by C. In cases where F is not 
small enough by itself, the system has to be sufficiently close to equilibrium in 
order to provide a valid gradient approximation through the smallness of the 
collision term C. Both the Kadanoff-Baym and the Botermans-Malfliet 
choice ( pT| ) are, of course, equivalent within the validity range of the first-order 
gradient approximation. Frequently, however, such equations are used beyond 
the limits of their validity as ad-hoc equations, and then the different versions 
may lead to different results. So far we have no physical condition to prefere 
one of the choices. The procedure, where in all Poisson brackets the Fin and 
Tout terms have consistently been replaced by /F and (1 =F /)F, respectively, is 
therefore optional. However, in doing so we gained some advantages. Beside 
the fact that quantum four-phase-space kinetic equation ( ^l| ) and the mass- 
shell equation are then exactly equivalent to each other, this set of equations has 
a particular features with respect to the definition of a nonequilibrium entropy 
flow in connection with the formulation of an exact H-theorem in certain cases. 
If we omit these substitutions, both these features would become approximate 
with deviations at the second-order gradient level. 

The equations so far presented, mostly with the Kadanoff-Baym choice 
( p^ ) , were the starting point for many derivations of extended Boltzmann and 
generalized kinetic equations, ever since these equations have been formulated 
in 1962. Most of those derivations use the equal-time reduction by integrating 

; the description to 
and refs. therein. 



6U27 2a,b9U3d 3]^,ti2 
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the four-phase-space equations over ener^^ff.^ 
three-phase-space information, cf. refs 
This can only consistently be done in the limit of small width F employing 
some kind of quasi-particle ansatz for the spectral function A{X,p). Partic- 
ular attention has been payed to the treatment of the time-derivative parts 
in the Poisson brackets, which in the four-phase-space formulation still ap- 
pear time-local, i.e. Markovian, while they lead to retardation effects in the 
equal-time reduction. Generalized quasiparticle ansatze were proposed, which 
essentially improve the quality and consistency of the approximation, provid- 
ing those extra terms to the naive Boltzmann equation (some times called 
additional collision term) which are responsible for the correct secoiiiLDrder 
virial corrections and the appropriate conservation of total energy, cf. ESo and 
refs. therein. However, all these derivations imply some information loss about 
the differential mass spectrum due to the inherent reduction to a 3-momentum 
representation of the distribution functions by some specific ansatz. With the 
aim to treat cases as those displayed in Figs. 2 and 3, where the differential 
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mass spectrum can be observed by di-lepton spectra, within a self-consistent 
non-equiblibrium approach, one has to treat the differential mass information 
dynamically, i.e. by means of Eq. ( po|) avoiding any kind of quasi-particle re- 
ductions and work with the full quantum four phase-space kinetic Eq. (|2|) . In 
the following we discuss the properties of this set of quantum kinetic equations. 

5 ^-derivable approximations 

The preceding considerations have shown that one needs a transport scheme 
adapted to broad resonances. Besides the conservation laws it should comply 
with requirements of unitarity and detailed balance. A practical suggestion has 
been given in ref.Q in terms of cross-sections. However, this picture is tied to the 
concept of asymptotic states and therefore not well suited for the general case, 
in particular, if more than one channel feeds into a broad resonance. Therefore, 
we suggest to revive the so-called ^-derivable scheme, originally proposed by 
BaymEj on the basia^of the generating functional, or partition sum, gijien by 
Luttinger and Wardtj, and later reformulated in terms of path-integralso. The 
auxiliary functional $ is given by two-particle irreducible vacuum diagrams. It 
solely depends on fully re-summed, i.e. self-consistently generated propagators 
iG{x,y) =< Tc'^{x)(p^ {y) >, where Tc indicates real-time contour ordering. 
The consistency is provided by the fact that <i> is the generating functional for 
the re-summed self-energy Yi{x,y) via functional variation of $ with respect 
to any contour ordered propagator G{y,x), i.e. 

-il](x,y) = T'5i$/<5iG(2/,a;). (23) 

An extension Id include classical fields or condensates into the scheme is pre- 
sented in ref. IIj In graphical terms, the variation (p3| ) with respect to G is 
realized by opening a propagator line in all diagrams of <&. The resulting set 
of thus opened diagrams must then be that of proper skeleton diagrams of E 
in terms of full propagators, i.e. void of any self-energy insertions. As a con- 
sequence, the ^-diagrams have to be two-particle irreducible, i.e. they cannot 
be decomposed into two pieces by cutting two propagator lines. 

The key property is that truncating the auxiliary functional $ to a limited 
subset of diagrams leads to a self-consistent, i.e closed, approximation scheme. 
Thereby the approximate forms of $ define effective theories, where $('*pp'' ) 
serves as a generating functional for the approximate self-energies 
through relation (^3|) , which then enter as driving terms Jor the Dyson's equa- 
tions of the different species in the system. As BaymliJ has shown, such a 
<I>-derivable approximation is conserving as related to global symmetries of the 
original theory. We explicitly cite the forms of the conserved Noether current 
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and of the energy-momentum tensor, cf. ref. 0, 



J 



(24) 

+ g^"' (^''^^ - £P°*) , (25) 



where 



-/:'"*(x) 



(5$ 



5X{X) 



is the density ofinteraction energy {\{X) locally scales the coupling strength of 
vertices, cf. ref.ll3) and the density of potential energy f p°* takes the following 
simple form within the first-order gradient approximation 



£P°*(X) - - 



ReE^ (ftf)+ ReG" (Fin t Tout) 



The first term of £p°* complies with quasi-particle expectations, namely mean 
potential times density, the second term displays the role of fluctuations / = 
Fin =p Tout in the potential energy density. This fluctuation term precisely 
arises form the Poisson bracket term in the kinetic Eq. ( |2l]) which induces a 
back-flow. It restores the Noether expressions ( p^ and ^25|) as being indeed 
the exactly conserved quantities. In this compensation we see the essential 
role of the fluctuation term in the quantum four-phase-space kinetic equation. 
Dropping or approximating this term would spoil the conservation laws. Before 
the gradient expansion, quantities ( p^ and (^5|) are exact integrals of equations 
of motion. While after the gradient expansion, they comply with the quantum 
four-phase-space kinetic equation ( |2l| ) up to the flrst-order gradient terms. 

At the same time the ^-derivable scheme provides thermodynamical con- 
sistency. The latter automatically implies correct detailed balance relations 
between the various transport processes. For multicomponent systems it leads 
to a actio = reactio principle. This implies that the properties of one species 
are not changed by the interaction with other species without affecting the 
properties of the latter ones, too. Some thermodynamic_fixamples have been 
considered recently, e.g., for the interacting ttTVA systemE3 and for a relativis- 
tic QED plasmaEj. 



6 Collision Term 

To further discuss the transport treatment we need an explicit form of the 
collision term (Hq) , which is provided from the $ functional in the — h matrix 
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notation via the variation rules ( p3| ) as 

C{X,p) = ^^F{X,p) - ^^FiX,p). (26) 
dF{X,p) oF{X,p) 

Here we assumed $ be transformed into the Wigner representation and all 
=FiG~+ and iG+^ to be replaced by the Wigner-densities F and F. Thus, 
the structure of the collision term can be inferred from the structure of the 
diagrams contributiiift to the functional To this end, in close analogy to the 
consideration of ref . we discuss various decompositions of the ^-functional, 
from which the in- and out-rates are derived. For the sake of physical trans- 
parency, we confine our treatment to the local case, where in the Wigner repre- 
sentation all the Green functions are taken at the same space-time coordinate 
X and all non-localities, i.e. derivative corrections, are disregarded. Deriva- 
tive corrections give rise to memory effects in the collision term, which will be 
analyzed separately for the specific case of the triangle diagram. 

Consider a given closed diagram of at this level specified by a certain 
number n\ of vertices and a certain contraction pattern. This fixes the topol- 
ogy of such a contour diagram. It leads to 2"^ different diagrams in the — h 
notation from the summation over all — h signs attached to each vertex. Any 
— h notation diagram of which contains vertices of either sign, can be de- 
composed into two pieces in such a way that each of the two sub-pieces contains 
vertices of only one type of sign^ 




with Fi • • • FmF{ • • • F4 linking the two amplitudes. The V*{X;pi, ...p[, ...) 
and V[3{X;pi^ ...p'l, ...) amplitudes represent multi-point vertex functions of 
only one sign for the vertices, i.e. they are either entirely time ordered (— 
vertices) or entirely anti-time ordered (-1- vertices). Here we used the fact that 
adjoint expressions are complex conjugate to each other. Each such vertex 
function is determined by normal Feynman diagram rules. Applying the matrix 

''To construct the decomposition, just deform a given mixed- vertex diagram of <3? in such a 
way that all + and — vertices are placed left and respectively right from a vertical division 
line and then cut along this line. 
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variation rules (g^) , we find that the considered <i> diagram gives the following 
contribution to the local part of the collision term (p6|) 

d^i dVi 



(27r)4 {2ttY 



(28) 



with the partial process rates 

RiX;p,,...p[,...)^ ^ Re{V:{X;pi,...p[,...)Vp{X;pi,...p[,...)}. (29) 

(a/3)e* 

The restriction to the real part arises, since with (a|/3) also the adjoint (/3|a) di- 
agram contributes to this collision term. However these rates are not necessar- 
ily positive. In this point, the generalized scheme differs from the conventional 
Boltzmann kinetics. 

An important example of approximate $ which we extensively use below 

is 



i$ 




(30) 



where logarithmic factors due to the special—features of the <i>-diagrammatic 
technique are written out explicitly, cf. ref. Ej. In this example we assume a 
system of fermions interacting via a two-body potential V — VoS{x — y), and, 
for the sake of simplicity, disregard its spin structure. The $ functional of Eq. 
( ^0|) results in the following local collision term 



C 



loc 



Pi d P2 



P3 



X 



I 



(27r)4 (27r)4 (27r)4 
xS^{p + pi-p2- Ps) {F2F3FF1 - F2F3FF1 



(31) 



where d is the spin (and maybe isospin) degeneracy factor. From this example 
one can see that the positive definiteness of transition rate is not evident. 

The first-order gradient corrections to the local collision term (28) are 
called memory corrections. Only diagrams of third and higher order in the 
number of vertices give rise to memory effects. In particular, only the last 
diagrajH of Eq. (^) gives rise to the memory correction, which is calculated 
in ref.l2il. 
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7 Entropy 



Compared to exact descriptions, which are time reversible, reduced description 
schemes in terms of relevant degrees of freedom have access only to some limited 
information and thus normally lead to irreversibility. In the Green's function 
formalism presented here the information loss arises from the truncation of the 
exact Martin-Schwinger hierarchy, where the exact one-uarticle Green function 
couples to the two-particle Green functions, cf. refs. cHa, which in turn are 
coupled to the three-particle level, etc. This truncation is achieved by the 
standard Wick decomposition, where all observables are expressed through one- 
particle propagators and therefore higher-order correlations are dropped. This 
step provides the Dyson's equation and the corresponding loss of information 
is expected to lead to a growth of entropy with time. 

We start with general manipulations which lead us to definition of the 
kinetic entropy flow. We multiply Eq. (|2l|) by -ln{F/A), Eq. (H) by 
(=F)ln(F/A), take their sum, integrate it over d^p/(27r)4, and finally sum the 
result over internal degrees of freedom like spin (Tr). Then wc arrive at the 
following relation 

d.^oM) = Tr I ^ In ^C{X,p), (32) 



where the quantity 

, dV A^r 

Tr ' ^ 



a{X,p) (33) 



(where (j{X,p) — =f[1 =f /] ln[l T /] — f^^f) obtained from the l.h.s. of the 
kinetic equation is interpreted as the local (Markovian) part of the entropy 
flow. Indeed, the Sy^^ has proper thermodynamic and quasiparticle limits O. 
However, to be sure that this is indeed the entropy flow we must prove the 
H-theorem for this quantity. 

First, let us consider the case, when memory corrections to the collision 
term are negligible. Then we can make use of the form (^Sj) of the local collision 
term. Thus, we arrive at the relation 

Tr / -^ln:^Goe(X,p)^Trl ( . . . . R 

J (27r)4 p ^°<=y 2 J (27r)4 (27r)4 

X {F,...Fi... h-Fi...} In^^^i27rf6^ (y^p, - Y^p'] . (34) 
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In case all rates R are non-negative, i.e. R > 0, this expression is non-negative, 
since (x — y) hi{x/y) > for any positive x and y. In particular, R > takes 
place for all $-functionals up to two vertices. Then the divergence of sj^^, is 
non-negative 

dMX) > 0, (35) 

which proves the _ff-theorem in this case with (^3|) as the nonequilibrium en- 
tropy flow. However, as has been mentioned above, we are unable to show that 
R always takes non-negative values for all <I>-functionals. 

If memory corrections are essential, the situation is even more involved. 
Let us consider this situation again at the example of the $ approximation 
given by Eq. ( ^0[ ) . We assume that the fermion-fermion potential interaction is 
such that the corresponding transition rate of the corresponding local collision 
term (^ij) is always non-negative, so that the iJ-theorem takes place in the 
local approximation, i.e. when we keep only (7'°*^. Here we will schematically 
describe calculations of ref.O which, to our opinion, illustrate a general strategy 
for the derivation of memory correction to the entropy, provided the i7-theorem 
holds for the local part. 

Now Eq. (|2|) takes the form 

d^sUX) - Tr / ^ m f + Tr / ^ In f C^, (36) 

where sj^^, is still the Markovian entropy flow defined by Eq. (^3|). Our aim 
here is to present the last term on the r.h.s. of Eq. (|3q ) in the form of full 
x-derivative 

^'Z (^ f = -5M-'^(i.cm(^) + Sc^^^X) (37) 

of some function s^^j^{X), which we then interpret as a non-Markovian cor- 
rection to the entropy flow of Eq. ( ^3| ) plus a correction (Sc^cm). For the 
memory induced by the triangle diagram of Eq.(pO[) detailed calculations of 
ref.Ej show that the smallness of the Scnicm, originating from small space-time 
gradients and small deviation from equilibrium, allows us to neglect this term 
as compared to the first term in r. h.s. of Eq. (I?!). Thus, we obtain 

d, {sL + 4\em) ^ Tr 1 ^ In > 0, (38) 

which is the //-theorem for the non-Markovian kinetic equation under consid- 
eration with s{^j, -I- s^^jom the proper entropy flow. The r.h.s. of Eq. (|3|) is 
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non-negative provided the corresponding transition rate in the local collision 
term of Eq. ( |3l| ) is non-negative. 

The explicit form of s^cm very complicated, see ref. Ej. In equilibrium 
at low temperatures we get sl^^^ - T3 InT which gives the leading correction 
to the standard Fermi-liquid entropy. This is the famous correction ESlO to 
the specific heat of liquid '^He. Since this correction is quite comparable (nu- 
merically) to the leading term in the specific heat T), one may claim that 
liquid ^ He is a liquid with quite strong memory effects from the point of view 
of kinetics. 

8 Summary 

A number of problems arising in different dynamical systems, e.g. in heavy-ion 
collisions, require an explicit treatment of dynamical evolution of particles with 
finite mass-width. This was demonstrated for the example of bremsstrahlung 
from a nuclear source, where the soft part of the spectrum can be reproduced 
only provided the mass-widths of nucleons in the source are taken explicitly 
into account. In this case the mass- width arises due to coUisional broadening 
of nucleons. Another example considered concerns propagation of broad res- 
onances (like /9-meson) in the medium. Decays of p-mesons are an important 
source of di-leptons radiated by excited nuclear matter. As shown, a consistent 
description of the invariant-mass spectrum of radiated di-leptons can be only 
achieved if one accounts for the in-medium modification of the p-meson width 
(more precisely, its spectral function). 

We have argued that the Kadanoff-Baym equation within the first-order 
gradient approximation, slightly modified to make the set of Dyson's equations 
exactly consistent (rather than up to the second-order gradient terms) , provide 
a proper frame for a quantum four-phase-space kinetic description that applies 
also to systems of unstable particles. This quantum four-phase-space kinetic 
equation proves to be charge and energy-^momentum conserving and thermo- 
dynamically consistent, provided it is based on a ^-derivable approximation. 
The $ functional also gives rise to a very natural representation of the collision 
term. Various self-consistent approximations are known since long time which 
do not explicitely use the <I>-derivable concept like self-consistent Born and 
T-matrix approximations. The advantage the $ functional method consists in 
offering a regular way of constructing various self-consistent approximations. 

We have also addressed the question whether a closed nonequilibrium sys- 
tem approaches the thermodynamic equilibrium during its evolution. We ob- 
tained a definite expression for a local (Markovian) entropy flow and were able 
to explicitly demonstrate the iJ-theorem for some of the common choices of $ 
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approximations. This expression holds beyond the quasiparticle picture and 

thus gcncraHzcs the well-known Boltzmann kinetic entropy. Memory effects in 
the quantum four-phase-space kinetics were discussed and a general strategy 
to deduce memory corrections to the entropy was outlined. 
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